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ABSTRACT 



Quantum stress-energy tensors of fields renormalized on a Schwarzschild background 
violate the classical energy conditions near the black hole. Nevertheless, the asso- 
ciated equilibrium thermodynamical entropy AS by which such fields augment the 
usual black hole entropy is found to be positive. More precisely, the derivative of AS 
with respect to radius, at fixed black hole mass, is found to vanish at the horizon for 
all regular renormalized stress-energy quantum tensors. For the cases of conformal 
scalar fields and U(l) gauge fields, the corresponding second derivative is positive, 
indicating that AS has a local minimum there. Explicit calculation shows that in- 
deed AS increases monotonically for increasing radius and is positive. (The same 
conclusions hold for a massless spin 1/2 field, but the accuracy of the stress-energy 
tensor we employ has not been confirmed, in contrast to the scalar and vector cases). 
None of these results would hold if the back-reaction of the radiation on the space- 
time geometry were ignored; consequently, one must regard AS* as arising from both 
the radiation fields and their effects on the gravitational field. The back-reaction, no 
matter how "small" , is therefore always significant in describing thermal properties 
of the spacetime geometries and fields near black holes. 
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I Introduction 

A black hole can exist in thermodynamical equilibrium provided that it is surrounded 
by radiation with a suitable distribution of stress-energy. In the semi-classical approach, 
such radiation is characterized by the expectation value of a stress-energy tensor obtained 
by renormalization of a quantum field on the classical spacetime geometry of a black hole. 
One can use such a stress-energy tensor as a source in the semi-classical Einstein equation, 

^ ^ renormalizedi (l) 

to calculate the change effected by the stress-energy tensor in the black hole's spacetime 
metric. This is the "back-reaction" problem associated with the spacetime geometry of a 
black hole in equilibrium. 

In this paper we use solutions of back-reaction problems of the above type to compute 
the thermodynamical entropy AS by which quantum fields augment the usual Bekenstein- 
Hawking black hole entropy Sbh = (1/4)A//^^^, where Ah is the area of the event horizon 
(Units are chosen such that G = c = kB = 1, but h ^ 1.). We consider explicitly the case of 
a Schwarzschild black hole surrounded by either a massless conformal scalar field or a U(l) 
gauge field (Maxwell field). (A massless spin 1/2 field is treated in the Appendix, but the 
accuracy of its stress-energy tensor has not to our knowledge been checked, in contrast to 
the conformal scalar and vector fields.) We show in all these cases that AS is positive. 

Our investigation shows rigorously that for all possible regular stress-energy tensors, the 
radial derivative of AS" vanishes at the horizon, for fixed black-hole mass; that is, AS" has 
there a local extremum with respect to radius. The form of the second derivative gives the 
criterion for a local minimum, which indeed occurs in all cases we have considered. Then by 
explicit calculation we show that AS is positive and monotonically increasing for increasing 
radius. Therefore the local minimum of AS* at the horizon is the only one and is its global 
minimum. As a consequence, the entropy is amenable to statistical interpretation. None 
of these features holds if the back-reaction of the fields on the spacetime metric is ignored. 
In this sense, AS must be regarded as arising from both the quantized radiation fields and 
from their effects on the gravitational field. 

We shall see, from the properties of the renormalized stress-energy tensors we employ and 
of the semi-classical Einstein equation, that we can obtain accurate fractional corrections 
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to the metric only in 0(e), where e = HM"'^, Mpi — h}^"^ is the Planck mass and M is the 
mass of the black hole. Because the usual black hole entropy Sbh — (47rM^)^~^ = 0(e~^), 
corrections to Sbh can be obtained in 0(e'') = 0(1) from fractional corrections of 0(e) in 
the metric. It turns out that these corrections are of the same order as the naive flat space 
radiation entropy (4/3)aT^V^, where a — (vr^/lS^^), Th — ^(SttM)"^ is the uncorrected 
Hawking temperature of a Schwarzschild black hole, and V is the flat space volume. Prom 
this fact alone it follows that the back-reaction cannot be ignored. 



II Stress-Energy Tensors 

Stress-energy tensors renormalized on a Schwarzschild background have been obtained 
in exact form for conformal scalar fields and for U(l) gauge fields, respectively, by Howard 
[1] and by Jensen and Ottewill [2]. Both results can be written in the form 

< Tjjf Prenormalized — P analytic ~l~ ^7j-2^^J^^4^ (^) 

where the analytic piece, in the case of a conformal scalar field, was given by Page [3]. The 
term is obtained from a numerical evaluation of a mode sum. The numerical piece is 
small compared to the analytic piece, and we do not include it in the calculations in this 
paper. This does not change any of our results qualitatively because both pieces separately 
obey the required regularity and consistency conditions. The analytic piece has the exact 
trace anomaly in both cases. 

The stress-energy tensors satisfy < >= on the Schwarzschild background with 
metric 



gi_,^ = diag 



2M. 2M , 2 2.2. 
^1 ) , r , r sm 



(3) 



r r 

These tensors represent the stress-energy distribution required to equilibrate the black hole 
with its own Hawking radiation. Each satisfies < >=< T^T > at the horizon r = 2M, 
which is required for regularity of the spacetime geometry [3] . Each has the asymptotic form 
of a fiat spacetime radiation stress-energy tensor at the uncorrected Hawking temperature 
at infinity of an ordinary Schwarzschild black hole, denoted here by Th — h{S7rM)~^. 

Dropping the angular brackets and displaying the analytic piece, one has for the confor- 
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mal scalar field [3] 

= -iaT^(^) (3 + 6w + 9w=^ + 12^2 + 15w^ + 18w^ - 99w^) , (4) 
r; = ^aT^(^) (1 + 2w + 3w^ + Aw^ + + 6w^ + 15w^) , (5) 

= = ^'^^■^^^^ (1 + 2W + 3^2 + 4«;3 + 5«;^ + - 9w^) , (6) 
where w = 2M/r. We have displayed the factor (1/2) explicitly because the scalar field has 
one helicity state while the vector field below has two. It is convenient in what follows to 
write 

where K = 38407r. For the U(l) vector field, we have [2] 

T* = -^aT^ (3 + 6W + 9w^ + 12w^ - 315w^ + 78w^ - 2A9w^) , (8) 

r; = ^aT^ (l + 2w + 3w^ - 76w^ + 295w^ - 54w^ + 285w^) , (9) 

= = ^aT^ (l + 2w + 3w^ + 44w^ - 305w^ + 66w^ - 579w^) . (10) 

In both cases > and the energy density — is negative in the vicinity of the event 
horizon, thus violating the weak energy condition. For the scalar field, the energy density 
is negative from r = 2M to r ~ 2.34M and for the vector field from r = 2M to r ~ 5.14M. 
Both tensors also violate the dominant energy condition in a region surrounding and bor- 
dering on the horizon. 



Ill Back-reaction on the Metric 

We obtain fractional corrections to the metric by setting 

9,u^gaAK + eK] (11) 

in the semi-classical Einstein equation (1). We work in linear order in e as required by 
V^T)^ = and V^((5G'(^) = 0, where SC^ is the Einstein operator linearized on a background 
satisfying = 0. The corrected geometry will be taken to be static and spherically sym- 
metric. Working out the equations as in [4], we find the corrected metric can be written 
as ^ 

ds^^-(l-'^I^]{l + 2ep{r))df+(l-'^^^ dP + r'du;', (12) 



where dcu'^ is the standard metric of a normal round unit sphere. To obtain m(r) and p{r) 
requires only simple radial integrals involving and T^. The angular components enter 
hnearized Einstein equations that hold automatically by virtue of V/^T^ = in a static 
spherical geometry. 

The mass function m(r) has the form 

m{r) = M{l + eij{r) + eCK-^), (13) 

with 

'"(^) = 47 r (-^/)47rf2c?f, (14) 
so ii{r) vanishes at the horizon. In (13), C is an undetermined integration constant that 
inspection of (12) shows is to be absorbed into M to obtain a renormahzed mass for the 
black hole. Thus, setting g'^'^ — shows that r = 2m = 2M(1 + eCK~^) — 2Mrenormaiized 
locates the event horizon. Note that, to the order we are working, we can write m(r) = 
M(l + eCK~^){l + e fi{r)) = Myeni}- + c/^(^))- The renormahzed mass will not be distin- 
guished notationally from the original Schwarzschild mass M in what follows, as the bare 
Schwarzschild mass has no physical meaning in the back-reaction problem. Therefore, we 
write 

m{r-) = M{l + eii{r)) = M + Mradir) (15) 

where, using (14), we see that Mj-^a — eM jiis the usual expression for the effective mass of 
a spherical source. 

For the scalar field, denoted where necessary by a subscript "s" , one finds [4] 

12 44 
Kiis = -{-w~^ + 2w-'^ + &w~^ -%\n{w) - 10w-6w^ + 22w^- — ). (16) 

2 3 3 

For the vector field, denoted by a subscript "v" , one finds [5] 

2 X 66 

Kn^ = -w-^ + 2w-2 + 6«;-^ - 81n(w) + 210^i; - 26^^ + - 248. (17) 

3 3 

In both (16) and (17), we note that the first term on the right, multiplied by eMK~^, gives 
the naive fiat-space value a T^V for radiation energy. 

The metric is completed by a determination of p which, like p, can be found from an 
elementary integration. Defining 

Kp = Kp + k, (18) 
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where A; is a constant of integration, we have 

p = - /" (T; - r/)(f - 2M)-^47rf2 dr. (19) 
e J2M 

For the scalar field, one finds [4] {Kps — Kpg + ks) 

1 /2 , 40 9 28 o 84\ 

i^Ps = 2 (3^" + - 8 - y w - lOw^ - — + — J . (20) 

Note that at the horizon r = 2M, or -u; = 1, we have Ps(l) = 0. The constant k for the 
scalar (vector) is denoted kg {ky) and will be determined below by a boundary condition. 
Similarly, for the vector field we have K p^j — K p^ + ky, where [5] 

2 40 

X p„ = -w-^ + 4w-^ - 8 ln(w) + —w + lOw^ + iw^ - 32, (21) 

3 3 

and Pt,(l) = at w = 1. 

Because both radiation stress-energy tensors are asymptotically constant, it is clear that 
the system composed of black hole plus equilibrium radiation must be put in a finite "box" . 
Otherwise, the fractional corrections e /i" to the metric would not remain small for sufficiently 
large radius. Physically, this means that the radiation in a box that is too large would 
collapse onto the black hole, producing a larger one. Hence, we must choose the radius Tq 
of the box such that it is less than the second positive root r* for r in g^^ — (the first 
zero corresponds to the horizon r = 2M). We shall also assume that the box radius Vo 
is sufficiently large that the stress-energy tensors we employ, which were constructed for 
infinite asymptotically flat spacetime, are a good approximation. Clearly, a flnitc radius 
would cut out some of the radial modes that were used in these calculations. However, if Vo 
is somewhat greater than the longest wavelength characteristic of Hawking radiation, which 
in turn is associated with the least-damped quasi- normal mode of lowest angular momentum 
for the field in question, then this effect should be neghgible. This wavelength A* is about 
42M for the conformal scalar field and is smaller for the higher-spin masslcss fields. Also, if 
To > A^, the explicit nature of the walls of the box {e.g., adiabatic versus diathermic) should 
not be important. For these reasons we shall assume throughout the remainder of this work 
that A* < To < r*. (Of course, one must also assume that M ~ Mpi, in any treatment 
based on (1).) If the radius rg were to approach the horizon, then explicit size and boundary 
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effects would fiave to be taken into account in tlie construction of < >, as shown in the 
work of Elster [6,7]. 

One convenient way to fix the constants kg and k^ is to impose a microcanonical boundary 
condition [4]. We fix To and imagine placing there an ideal massless perfectly reflecting wall. 
Outside To, we then have an ordinary Schwarzschild spacetime 

ds^^-(l- dt^+(l- ~\r^ + rV, (22) 



for r > Tq. Continuity of the three-metric induced by metrics (12) and (22) on the world 
tube r — To fixes the constant /c, i.e., kg or /c^, in p by the relation 

k = -Kp{ro). (23) 

There are finite discontinuities in the extrinsic curvature of the world tube r = Tq [4], 
but these, and other properties of the box wall, are of no interest in the present analysis, 
as we argued above. The spacetime geometry, including back-reaction, is now completely 
determined by (22) for r > Tq, and for r < Tq by 

ds^^-U- [1 + 2e (p(r) - p{To))]de + f 1 - + du\ (24) 
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IV Temperature 

If we release a small packet of energy from a closed box containing a black hole through 
a long thin radial tube, it will undergo a red-shift and approach the asymptotic temperature 



T = 



27r ' 



(25) 



where kh is the surface gravity of the event horizon. For an ordinary Schwarzschild black 
hole (ignoring the radiation), one finds kh — (4M)~^ and Too = Th = h{87rM)~^. However, 
the stress-energy of the radiation changes the surface gravity of the horizon to 

1 



AM 



1 + e(p - /i) + 87rr^ Tl 



\r=2Mi 



(26) 



as a straightforward calculation shows [4]. With the microcanonical boundary conditions, 

we can use (23) to obtain from (25) and (26) 

where n takes the value Ug — 12 for the scalar field and riy — 304 for the vector field. The 
local temperature at the boundary of the box is obtained by blue-shifting (27) from infinity 
back to To- We find from 

Tioc^T^[-gu{ro)V^ (28) 

that 



Tioc{ro) 



SttM 



1 — e p(ro) -l- eriK ^ 



1 - 



2m(r„) 

To 



n-1/2 



(29) 



The temperature TJoo unlike Too is actually independent of the boundary condition that 
determines the constant k, as explained in detail in [4]. Indeed, it can be readily verified 
by the reader that k cancels out in 0(e) in the expression (28) for Tioc- Either measure of 
temperature, T^o or Tioc, can be used to calculate the same entropy in conjunction with an 
appropriate measure of energy. This is quite important: it means that the specific boundary 
condition chosen does not affect the calculated entropy, as we shall see below. 
V Thermodynamical Entropy 

One way to calculate the entropy is as follows. Fix the radius of a closed box. The 
measure of energy in the box conjugate to the asymptotic inverse temperature Poo = 
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is then the Arnowitt-Deser-Misner (ADM) mass m(ro) determined at spatial infinity. The 
first law of thermodynamics for slightly differing equilibrium configurations tells us that 



dS — /3oo dm {dvo — 0), 



(30) 



where S{ro) is the total entropy in the box. By this method we seem to obtain only the 
total entropy S{ro) rather than the distribution of entropy in the given box, <S'(r), for r < Tq, 
where S{r) denotes the total entropy inside the radius r. However, the latter can be obtained 
by using the quasi-local energy E [8-11], which for static spherical metrics like those treated 
here is given for any radius r < by 



E{r) = r- r[/^(r)]^/^ 



(31) 



with g^^{r) determined by (24), the metric for r < To- This energy, unUke m, does not depend 
on asymptotic flatness in its definition, nor even on the existence of an asymptotically flat 
region [10,11]. Furthermore, even the "normalization" of the zero of energy [10,11] that is 
incorporated in E as given in (31) does not affect the calculated entropy, as it certainly 
should not. (This "normalization" is intended to make E approach the ADM mass in an 
asymptotically flat region, if such a region exists.) Similarly, the inverse local temperature 
/?(r) = Ti~^{r), r < ro, is independent of the boundary conditions as mentioned above. 
Hence, the value of the entropy depends neither on the zero of energy nor on the existence 
of an asymptotic region. 

Therefore, to obtain S{r), in place of (30) we can write 



dS^PdE {dr^O,r<ro). 



(32) 



Choosing M and r as independent variables, and fixing r, we can readily integrate (32) to 
obtain 5" up to a function of r and a constant. Prom (29) we have 

2m(r)"' 



StiM 

h 



1 + e p{r) — euK 



-1 



(33) 



and from (15), (24), and (31), holding r fixed. 



dE 



1 -e/i + eM 



dM 



- -1/2 


^ 2m(r) 








r 



n-l/2 



dM. 



(34) 
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One can see directly for any r < To that (5oodm — (5dE where, of course, one replaces 
by r in the formulas for (5oo and m to establish this result. This equality means that we 
can calculate S{r) for any r < Tq. The key point of this discussion is that one can think of 
adding layer upon layer of entropy, associated with the black hole and a given <Tl^ > that 
is valid from r = 2M io r — ro, beginning at r = 2M and ending at r = Tq. (Additivity of 
entropy in configurations analogous to this case is established in [12], but our method here 
establishes it independently.) 

Observe that from fractional changes of 0(e) in the metric, which affect the surface 
gravity and temperature in this order, we are able to calculate from (32) departures of 
0(e°) = 0(1) from the usual black hole entropy Sbh = (47rM^)^~^ = Ane~^. But in fact all 
of the corrections to the entropy are of the same order as the naive flat-space entropy itself: 

Lt'hV = ^(^)(7r^)'(^7rr3) = ^(^w'^) - 0(1) x w-\ (35) 

The in (35) cancel out, leaving only a function of w = 2Mr~^. 
Combining (33) and (34) yields 



dS = ^dM + Stt 
n 



aw 



dw, (36) 



with dr — 0. Integration of (36) gives an expression of the form 



S = — ^ + /:^S{w) + /(^), {l<w<Wo^ 2M/ro) (37) 

where the first term is the usual Bekenstein-Hawking expression Sbh for the black hole 
entropy, the second term is a function of w determined up to an additive integration constant 
by the second term on the right of (36), and / is a dimensionless function of r that does not 
depend on M. The appearance of a function / in (37) can be understood as follows. Since 

1 /2 

our problem involves three mass or length scales Mpianck = ^ j the mass of the black hole, 
M, and a radius r < rg, there are, for a given r, exactly three dimensionless parameters 
one can define, namely, e = HM"'^, w — 2M/r and r/h}^^. However, the first two terms on 
the right of (37) depend only on e and w, respectively. Thus, if the entropy S depends on 
r/fi^^'^, it can only do so through a separate function of this parameter. 

Let us first dispose of the dimensionless function /, which clearly can depend only on 
{r/h^^"^), where h^^"^ is the Planck length in our units. It seems that such a term could only 
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arise in a theory taking quantum gravity into explicit account because the semi-classical 
theory has incorporated the dimensionless terms involving h/M"^ and 2M/r. (Of course, 
quantum gravity could modify terms of these latter two types quantitatively) On dimen- 
sional grounds, therefore, we take / = in the semi-classical theory. (A formal argument 
that f — based on [9] can be constructed [13].) The possibility of an additive constant 
will be discussed when we treat AS below. 

In considering AS, which will be given explicitly below, we first note the significant 
property that 



d{AS) 
dw 



Stt 



-1/ \ T^-l -1 

w (P — U) + nK w 

aw 



(38) 



vanishes at the horizon w — 1. Therefore, for a fixed black hole mass M, the derivative with 
respect to r of AS vanishes at the horizon. Thus AS has a local extremum with respect 
to r at the horizon. This result follows from several general features that will be enjoyed 
by all regular renormalized stress-energy tensors on the Schwarzschild background and the 
back-reactions they induce, not just the cases analyzed here. First, /x vanishes at the horizon 
by virtue of the black hole's mass having been suitably renormalized. Second, p vanishes at 
the horizon, as follows from (19) and the regularity condition T/ = T^** at the horizon [3]. 
More precisely, we have that 



lim„ 



r*-r: 



,1+ [-^ — exists. (39) 

Third, the last two terms on the right of (38) add to zero at the horizon because there 
the Hamiltonian constraint {G\ — SnT^ — 0) holds. Furthermore, note that if the fractional 
effects of 0(e) in the temperature induced by the back- reaction were neglected, the derivative 
(38) would not vanish at the horizon, a property that the reader can verify. 

Is the local extremum of A^* at the horizon a local minimum? To answer this we calculate 
d\AS) 



' dw dw dw"^ 



which becomes, at the horizon w = 1, 

&\AS)^ 
dw"^ 

or , equivalently, with M fixed. 



kt)=l 



\ dw dw'^ 



\w=l 



d\A S) 



\r=2M 



327r^M^ 

n 



AM 



dr 



8t; - ( 



1 - 2M/r 



\r=2M- 



(40) 



(41) 



(42) 
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Hence we need only examine the stress-tensors. In all the cases we consider (conformal scalar, 
vector, massless fermion), (41) and (42) are positive so that A^" takes a local minimum with 
respect to radius at the horizon. This suggests, but does not prove, that AS is non-negative. 

The local minimum of AS* at the horizon and the fact that Sbh in the expression (37) 
for the total entropy S contains the renormalized mass M of the hole motivate the choice of 
the remaining additive constant in AS, which can only be a pure number, to be such that 
AS — a,t w — 1. For w — 1, with no "room" for the fields to contribute anything further, 
one then obtains only the Bekenstein-Hawking entropy (l/4)Ai/^~^, as would be expected. 
With the choice AS{w = 1) = 0, we obtain for the conformal scalar field [14,15] 

87r,l,/8 , 8 2 o 1 32, , , 40 „ o 104 , 16\ 



" 2^ U"^" ^ 3^^ ^ ^ y ^""^"^^ ~ y ~ + t ~ yj ^^^^ 

for 1 > ty > Wo- Similarly, for the electromagnetic or U(l) gauge field we find 

Stt /8 o 8 2 o 1 / N 40 „ 2 344 496\ ,,,, 

AS^ = — -w-^ + + 8w-^ - 96 ln(w) H w - 8w^ H . (44) 

K \9 3 ^ ^ 3 9 9 / ^ ^ 

In both expressions, the naive flat-space radiation entropy term (35) appears as the first 
term on the right. Both ASg and A^*^ are positive for 1 > w > > = 2Mr~^ and 
vanish at w — 1. Hence, in that they are positive, both are amenable to arguments relating 
thermodynamical and statistical entropy. It has not heretofore been evident that this desir- 
able feature would be present in the semi-classical theory. The reader can verify, by omitting 
the back-reaction terms in the inverse temperature (33), that not only is the vanishing slope 
of AS" at w = 1 lost, but also that the value of the resulting "AS"' , normalized as above, is 
no longer positive for the range 1 > w > Wo- In this fundamental sense, we conclude that 
the back-reaction, however small quantitatively in its effects on the metric near a black hole, 
can never be regarded as negligible. 
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Appendix 

Here we outline the calculation of A^" for a massless spin 1/2 field. We use the stress- 
energy tensor given in [16]. As far as we have been able to determine, its accuracy has not 
been verified by an exact numerical analysis, unlike the two cases wc treated in the body of 
the text. This tensor has also been used in a calculation similar to the one presented here 
in [17], where qualitatively different results were obtained for the entropy AS. 

The stress-energy tensor is given by 

= -laT^ il) (3 + 6W + 9w^ + 12w' + ^^w^ + - 69w'') , (^1) 

3 8V 77 / 

= il) {l + 2w + 3w' - - - + fw') , (.12) 

= r; = Ian il) {l + 2w + - fw' - fw' + f + 23^^) . (.13) 
We find for n and p 

7 /2 . 2 1 , / X 90 62 2 46 , 16\ , , 
Kpf = - -w~-^ + 2w~^ + 6w~^ - 81n(w) - —w - —w^ + —W^ - 17 > (^4) 

8 V3 7 7 3 7 / 

7 /2 _2 . -1 ow N 200 50 2 52 3 32\ 
Kpf = 8 (3^ + 4^^^ - 81n(w) - —w - —w^ + —w^ + yj , (^5) 

where the subscript "/" denotes "fermion". The formulas for temperature and inverse 
temperature have the same form as before with Uf = —4. The quantity /S.S enjoys all the 
same basic properties as for the conformal scalar and vector fields. It is given by 



87r,7, /8 o 8 2 . 1 24, , , 200 56 2 800 o 424 



63 63 



(A6) 

and is positive. 
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